QED Model of Resonance Phenomena in Multichain Qubit Systems 
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The analytical solution of the task of the interaction of quantized EM-field with multichain qubit 
system by taking into account both the intrachain and interchain qubit coupling has been obtained 
for the firsr time. It is argued the applicability of SSH-model for perfect carbon zigzag nanotubes. 
The appearance in optical spectra of additional lines, which are consequence of quantum nature of 
EM-field, leading to Rabi wave packets formation, is predicted. 
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Quantum electrodynamics (QED) consideration of in- 
teraction of electromagnetic (EM) field with matter is 
used at present very rare in practical applications, for 
instance, quantum nature of EM-field is not taken into 
account by its interaction with nanoobjects like to car- 
bon nanotubes, both in theoretical and experimental as- 
pects Q, Qj. There exist at present in QED-theory the 
analytical solutions of the task of the interaction of quan- 
tized EM-field with one qubit, Jaynes-Cummings model 
(JCM), with multiqubit systems without interaction be- 
tween qubits, Tavis-Cummings model, and the model 
described in |3j, generalizing Tavis-Cummings model by 
taking into account the ID-coupling between qubits. The 
model to be considered continues given series, that is, it 
is the model of the interaction of quantized EM-field with 
multichain qubit system by taking into account both the 
intrachain and interchain qubit coupling. Having in view 
the practical applications, the zigzag-shaped carbon nan- 
otubes (ZSCNTs) will represent multichain qubit system 
in the task aimed. 

ZSCNTs can be considered to be the set of carbon 
backbones of trans-polyacetylene (t-PA) chains, which 
are connected between themselves. It is substantial, that 
the theory of electronic properties of t-PA chains is very 
well developed and it has very good experimental con- 
firmation. Let us touch on given subject in more detail. 
Pioneering works [4j, [5| were opening new era in the 
physics of conjugated ID-conductors. The authors have 
found the most simple way to describe mathematically 
rather complicated system - the chain of t-PA, vibronic 
system of which represents itself the example of Fermi liq- 
uid. It is now well known SSH-model. The most substan- 
tial suggestion in SSH-model is the suggestion, that the 
only dimerization coordinate u n of the n-th CH -group, 
n= 1, N along chain molecular-symmetry axis x is sub- 
stantial for determination of main physical properties of 
the material. Other five degrees of freedom were not 
taken into consideration. Nevertheless, the model has 
obtained magnificent experimental confirmation. Natu- 
rally, it is reasonable to suggest, that given success is 
the consequence of some general principle. Given gen- 
eral principle is really exists and main idea was proposed 
by Slater at the earliest stage of quantum physics era 



already in 1924. It is - "Any atom may in fact be sup- 
posed to communicate with other atoms all the time it 
is in stationary state, by means of virtual field of radia- 
tion, originating from oscillators having the frequencies 
of possible quantum transitions..." |6]. Given idea ob- 
tains its development, if to clarify the origin of virtual 
field of radiation. It is shown in 7], that Coulomb field 
in ID-systems has the character of radiation field and it 
can exist without the sources, which have created given 
field. Consequently, Coulomb field can be considered to 
be "virtual" field in Slater principle. It produces pref- 
erential direction in atom communication the only along 
chain axis, and it explains qualitatively the success of 
SSH-model. 

It is remarkable, that SSH-model along with the phys- 
ical basis of the existence of solitons, polarons, breathers, 
formed in 7r-electronic subsystem (7r-solitons, 7r-polarons, 
7r-breathers) contains in implicit form also the basis for 
the existence of similar quasiparticles in cr-electronic sub- 
system, that is cr-solitons, er-polarons, er-breathers, dis- 
cussed in [8J. SSH-CT- polarons have been experimentally 
detected in carbynes [8[ , where the formation of polaron 
lattice was established. Slater principle allows to sug- 
gest, that the picture in perfect zigzag NTs can be quite 
similar, that is, despite the strong interaction between 
the chains, the physics of propagation processes will be 
very similar to that in t-PA and in carbynes. Given sug- 
gestion has experimental confirmation, it follows from 
ESR-studies in rather perfect CZSNTs, produced by high 
energy ion modification (HEIM) of diamond single crysa- 
tals @. In other words, Peierls transition is taking place, 
that is, perfect zigzag NTs , characterized by (m, 0) in- 
dices, will have bandgap like to semiconductors at any 
m, including the case m = 3q, q £ N, for which the 
theory, existing at present, predicts the metallic proper- 
ties. Given difference is consequence of that, all known 
to our knowledge theoretical works do not take Peierls 
transition into account, see for instance [1], 0]. To ob- 
serve quantum coherent effects on NTs, including Peierls 
transition, the ensemble of NTs has to be homogeneous. 
Moreover, to satisfy a requirement of onedimensionality 
and homogeneity any dispersion in axis direction, chiral- 
ity, length and especially in diameter both for single NT 



along its axis and between different NTs in ensemble has 
to be absent, axial symmetry has to be also retained, that 
is, there are additional requirements in comparison with, 
for examplle, t-PA technology. The technology, based on 
HEIM, discussed in review [9], satisfy given requirements. 
The CVD-technology and many similar ones seem to be 
not satisfying to above-listed requirements at present. It 
means, that experimental results and their theoretical 
treatment will be different in both the cases, that really 
takes place, if to compare ESR results. 

We can now consider the task put by. In the frames of 
SSH-model ZSCNT represents itself autonomous dynam- 
ical system with discrete circular symmetry consisting of 
finite number n € N of carbon backbones of t-PA chains, 
which are placed periodically along transverse angle co- 
ordinate. Longitudinal axes {xi},i = l,n, of individual 
chains can be directed both along element of cylinder and 
along generatrix of any other smooth figure with axial 
symmetry. It is taken into account, that in the frames of 
SSH-model, the adjacent chains, which represent them- 
selves a mirror of each other, will be indistinguishable. 
Here we will consider the only SSH-er-polarons to be opti- 
cally active centers, which produce polaron lattice. Each 
(T-polaron in accordance with experiment [8j can be ap- 
proximated like to guantum dot in |3j by two-level gu- 
bit. Then the Hamiltonian, proposed in given work can 
be generalized. The insufficient for the model local field 
term was omitted. (Local field term seems to be playing 
minor role by description of cr-polarons in comparison 
with quantum dots, since size of quantum dots is greatly 
exceeding the size of cr-polarons). We will use the ap- 
paratus of hypercomplex n-numbers. Let us remember, 
that hypercomplex n-numbers are defined to be elements 
of commutative ring 
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that is, it is direct sum of n fields of complex numbers C, 
n € N. Then the QED-Hamiltonian, considered to be hy- 
percomplex operator n-number, for cr-polaron system of 
zigzag NT, consisting of n t-PA chains and connected be- 
tween themselves in that way, in order to produce rolled 
up graphene sheet, and interacting with EM-field in ma- 
trix representation is 



[H] = \H a \ + [H F ] + \H aF \ + [TiaaY 



(2) 



The rotating wave approximation and the single-mode 
approximation of EM-field are used. All the components 
in ^ are considered to be hypercomplex operator n- 
numbers and they are the following. [Ho] represents the 
operator of cr-polaron subsystem energy in the absence 
of interaction between cr-polarons themselves and with 
EM-field. It is 



[n,] = {Hw /2)J2T,^A^y 
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tr-polaron in j-th chain. The second term 
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[H F ] = ^ E « + aM J 
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is the Hamiltonian of the free EM-field, which is repre- 
sented in the form of hypercomplex operator n-number. 
The component of the Hamiltonian ([2]) 
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describes the interaction of cr-polaron sybsystem with 
EM-field, where g is the interaction constant. The Hamil- 
tonian 
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where hR' !, are the energies, characterizing intrachain 
(I = j) and interchain (I ^ j) polaron-polaron interac- 
tion for the excited (£^) and ground (^ 2 -*) states of j-th 
chain, \b m j) , \cimj) are ground and excited states corre- 
spondingly of m-th cr-polaron of j-th chain. In ^ to © 
[eip is j-th power of the circulant matrix [ei], which is 
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Hamiltonian in the form like to (JBJ) at n = 1 is usually 
used for description of tunneling between the states with 
equal energies, inparticular, for tunneling between quan- 
tum dot states [3|. Hamiltonian © at any n describes 
actually the connection between pairs of the states, which 
satisfy the following condition - the first state in any pair 
results from the second state (and vice versa) by time 
reversal. It is known, that for given states Cooper effect 
takes place. Therefore, the application of Hamiltonian 
like to ([HI) is possible for any pair of time reversal sym- 
metric states with equal energy. By the way, if one omits 
last term in Hamiltonian ((2), it goes into n-chain gen- 
eralization of well-known Tavis-Cummings Hamiltonian. 
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Figure 1: Temporal dependence of the integral inversion in 
the QD chain for a coherent initial state of light ((n) — 4). 
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Figure 2: Spectral dependence of the signal amplitude, corre- 
sponding to temporal dependence, presented in Fig.l 



We can represent the state vector of the "NT+EM- 
field" system in terms of the eigenstates of isolated po- 
laron and photon number states in the form 
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Here, |6 mi ,/) = |6 mj -) ® |i), |o m j,/) = \a m j) ® |0> where 
|Z) is the EM-field Fock state with I photons, A 3 ml (t), 
B J m ,(t) are the unknown probability amplitudes. Let 
us pay the attention, that [|\I/(£))] is one state and the 



representation in the form of sum over j is like to repre- 
sentation of complex- valued wave function ^(t) in the 
form W(t) = ^i(t) + i^ 2 {t) of two real- valued func- 
tions. Then, solving the nonstationary hypercomplex 
Schrodinger equation for the state vector [|\l/ (£))], we ob- 
tain in continuum limit 
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[$ (h,0)]exp{it([9 l (h)] - gy/T+T\x])}e ihx dh 
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where x is hypercomplex axis x = [x, x, ..., x], [$ (a;, £)] is 

[*W)J =e^^^M exp ^ [ ^ (a; , i)]] (10) 

in which [er z ] is Pauli z-matrix, A = r _1 is phenomenolog- 
ical relaxation factor. In its turn [\l/'(a;,i)] is continuous 
limit of functional block matrix of discrete variable m, 
which is 
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[A m ,i{t)] 
[B m ,i+i(t)] 
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that is, it is consisting of two [n x n] matrices of proba- 
bility amplitudes 



3=0 

ra-1 

[£m,Z+l(*)]=E<*+l(*)N J 
J=0 



(12) 



which are determined by relationship ([5]). Further, ma- 
trix [9(h)] in © is 



{o(h)] = ^{([e 1 (h)] + [e 2 (h)])0[E 2 }} 
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where [E 2 ] is unit [2 x 2]-matrix, [61(h)] and [6 I 2(^)] are 

(14) 
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[e 2 (h)} = ^ 2 }{2-a 2 (h-^) 2 } 

Here [£1], [£2] are [n x h] matrices of coefficients in ([6]), 
that is 
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where one takes into account, that in view of axial sym- 
metry [£j 2 ] do not depend on j. Consequently, we have 

[&] = Kii, U] = [6]- Matrix [ X ] in ® is 



M = [-^n] ® [o"x] exp-i[cr z ](o; - w )i. 
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Matrix elements of [$ l (x, t)] are 
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exp {i y j exp {$j{h) — gVl — lKj(h))}dh, 

3=0 



where 8' (h, Q),$j(h), Kj(h) are determined by eigenval- 
ues k Q € C, a = 0,n-l of $ z (/i,0), 0(/i) and%(/i), which 
are considered to be n-numbers. They are 
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Then the hypercomplex solution can be represented in 
the form of sum of n solutions for n chains, that is, hy- 
percomplex n- number &(x,t) is 



n-l 
$ / (^,t) = ^^( a;! t), 

where the solution for q-th chain <t l (x,t) is 



p=0 



(21) 
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in which the matrix elements & (x,t) are determined 
by (jTTJ) . The relationship (j2T)f by taking into account © 



- (fT9|) determines Rabi-wave packet, which propagates 
along individual chain of zigzag NT. It is also clear, that 
now the analysis of Rabi-wave packet dynamics for indi- 
vidual NT-component will be the same, that in [3|. Let 
us reproduce Fig. 4 from [3|, where temporal dependence 
of the integral inversion in the QD chain for a coherent 
initial state of light ((n) = 4) is presented. It is seen, 
that an originally Gaussian packet temporally oscillates, 
at that oscillations collapse to zero quickly, but revive 
with time increasing in, and what is characteristic, in 
another area of space j3| ■ 

Temporal dependence (Fig.l) of the integral inversion 
gives in implicit form the way for comparison of theoret- 
ical results with any stationary optical experiments. Re- 
ally, it is sufficient to make a Fourier transform of given 
temporal dependence (Fig. 2). It will be proportional 
to signal amplitudes of infrared (IR) absorbance, IR- 
transmittance, IR-reflectance or Raman scattering, since 
they arc determined by population difference. There- 
fore, the QED-model presented predicts the appearance 
of additional lines, for which Rabi waves are responsi- 
ble, in stationary optical spectra. It is understandable, 
that Rabi waves can be registered in any ID-systems by 
stationary optical measurements, if the electron-photon 
coupling is rather strong. The comparison with exper- 
iment is given in [lOJ , which confirms given theoretical 
results. 

Thus, in fact, we have obtained for the first time in 
QED the analytical solution of the task of the interac- 
tion of quantized EM-field with multichain qubit system 
by taking into account both the intrachain and inter- 
chain qubit coupling. The model considered continues 
the QED-model series, in which the next model general- 
izes previous model and which at the same time can be 
solved analytically, that is, it ranks on a par with Jaynes- 
Cummings model for the one qubit, Tavis-Cummings 
model, which is JCM-generalization for multiqubit sys- 
tems and the model described in 3], generalizing Tavis- 
Cummings model by taking into account the ID-coupling 
between qubits. 
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QED-model for the multichain qubit system with interactions of qubits and chains between them- 
selves on the example of the system of a-polarons in carbon zigzag nanotubes, interacting with 
quantized EM-field, is considered analytically. The possibility of experimental detection of Rabi 
waves in conventional stationary optical experiments for any quasi- ID system with strong electron- 
photon interaction is predicted. 
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Quantum electrodynamics (QED) consideration of in- 
teraction of electromagnetic (EM) field with matter is 
substantial for many practical applications, for instance 
for implementing of quantum computation, and it is 
studied extensively in two main directions. They are 
the Jaynes-Cummings model (JCM) extending by con- 
sidering of multiple atoms and multiple mode field and 
JCM-extending by introducing of nonlinearity in JCM- 
Hamiltonian. The examples of the synthesis of given 
two directions are recent works [1], 0], where Tavis- 
Cummings model [3| is developed by taking into account 
the interaction between qubits on the one hand and si- 
multaneously a nonlinearity is also included into con- 
sideration on the other hand. The authors predict the 
new coherent effect of nonlinear quantum optics - spa- 
tial propagation of Rabi oscillations (Rabi waves) in ID 
quantum dot chain. It represents the interest in a QED- 
consideration of the interaction with EM-field of multi- 
chain systems by including the interchain interaction. It 
is aim of the paper presented. Especially interesting to 
consider practically significant case of carbon nanotubes 
(NT) to be a matter subsystem of the system matter 
+ quantized EM-field. Really, for instance, NT of zigzag 
kind and NT of armchair kind can be considered to be the 
set of trans-polyacetylene (t-PA) and cis-polyacetylcnc 
(c-PA) chains correspondingly, which are connected be- 
tween themselves. There seems to be the most interesting 
to consider zigzag carbon NT, since the theory of elec- 
tronic properties of t-PA chains is very well developed 
and it has very good experimental confirmation. Let us 
touch on given subject in more detail. Pioneering works 
[J| [5( were opening new era in the physics of conjugated 
ID conductors. The authors have found the most sim- 
ple way to describe mathematically rather complicated 
system - the chain of t-PA, electronic system of which 
represents itself the example of Fermi liquid. It is now 
well known SSH-model. The most substantial suggestion 
in SSH-model is the suggestion, that the only dimeriza- 
tion coordinate u n of the ra-th C-ff-group, n = l,N along 
chain molecular-symmetry axis x is substantial for de- 
termination of main physical properties of the material. 
Other five degrees of freedom were not taken into con- 



sideration. Nevertheless, the model has obtained mag- 
nificent experimental confirmation. Naturally, it is rea- 
sonable to suggest, that given success is the consequence 
of some general princip. Given general princip is really 
exists and main idea was proposed by Slater at the ear- 
liest stage of quantum physics era already in 1924. It is 
- "Any atom may in fact be supposed to communicate 
with other atoms all the time it is in stationary state, by 
means of virtual field of radiation, originating from oscil- 
lators having the frequencies of possible quantum tran- 
sitions..." [6(. Given idea obtains its development, if to 
clarify the origin of virtual field of radiation. It is shown 
in [7], that Coulomb field in ID systems has the charac- 
ter of radiation field and it can exist without the sources, 
which have created given field. It can be considered to be 
virtual field in Slater princip. It produces preferential di- 
rection in atom communication the only along chain axis, 
and it explains qualitatively the success of SSH-model. 
It is remarkable, that SSH-model contains in implicit 
form along with the physical basis of the existence soli- 
tons, polarons, breathers, formed in 7r-electronic subsys- 
tem (7r-solitons, 7r-polarons, 7r-breathers), also the basis 
for the existence of similar quasiparticles in er-electronic 
subsystem, that is cr-solitons, er-polarons, cr-breathcrs. 
The cause is the same two-fold degeneration of ground 
state of the whole electronic system, energy of which 
has the form of Coleman- Weinberg potential with two 
minima at the values of dimerization coordinate uq and 
— uq. SSH-cr-polarons have recently been experimentally 
detected in related material - carbynes [9] . where the for- 
mation of polaron lattice was proposed. Corresponding 
chain state is optically active and it is characterized by 
the set of lines in IR-spectra, which were assigned with 
new phenomenon - antiferroelectric spin wave resonance 
(AFESWR). Slater princip allows to suggest, that the 
picture in zigzag NT can be quite similar, that is, despite 
the strong interaction between the chains, the physics of 
propagation processes will be very similar to that in t- 
PA and carbynes and propagation process in zigzag NT 
will represent the sum of n individual processes. We can 
now return to the problem put by. In the frames of SSH- 
model zigzag NT represents itself autonomous dynamical 



system with discrete circular symmetry consisting of fi- 
nite number n £ N of t-PA chains, which are placed peri- 
odically along transverse angle coordinate. Longitudinal 
axes {xi},i — l,n, of individual chains can be directed 
both along element of cylinder and along generatrix of 
any other smooth figure with axial symmetry. It is taken 
into account, that in the frames of SSH-model, which 
take into consideration the dynamical processes the only 
along chain direction, the adjacent chains, which repre- 
sent themselves a mirror of each other, will be indistin- 
guishable. Here we will consider the only SSH-cr-polarons 
to be optically active centers, which produce polaron lat- 
tice. Each cr-polaron in accordance with experiment [9j 
can be approximated like to guantum dot in [1] by two- 
level gubit. Then for any individual chain the Hamilto- 
nian, proposed in [2J can be used. Further to take into 
consideration all the rest chains and the interaction be- 
tween them, we will use the apparatus of hypercomplex 
n- numbers. Let us remember briefly, that hypercomplex 
n-numbers are defined to be elements of commutative 
ring 



Z n = C(BC< 



©C, 



(1) 



that is, it is direct sum of n fields of complex numbers 
C, n £ N . It means, that any hypercomplex n-number 
z £ Z n is n-dimensional quantity with the components 
k a £ C, a = 0, n — 1, and it in row matrix form z £ Z n 



z = [k a ] = [ko,ki,...,k n -i] 
It can be represented in the form 



n-1 

E 

a=0 



^a^a. ) 



(2) 



(3) 



where 7r Q are basis elements of Z n (and simultaneously 
basis elements of the linear space of n-dimensional lines 
and n-dimensional row matrix). They are 



(4) 



7ro = [l,0,...,0,0],7r 1 = [0,l,...,0,0], 

...,7T„_l = [0,0,...,0,l]. 

Basis elements ix a possess by projection properties 

n-1 

TT a TT a = lT a 8 a f3, 2_^ TT a = 1, ZTT a — k a lT a (5) 

a=0 



In other words the set of k a £ C,a = 0, n — 1 is the set 
of eigenvalues of hypercomplex n-number z £ Z n , the set 
of {7r a }, a — 0,n — 1 is eigenbasis of Z n -algebra. Then 
the QED-Hamiltonian, considered to be hypercomplex 
operator n-number, for cr-polaron system of zigzag NT, 
consisting of n t-PA chains and connected between them- 
selves in that way, in order to produce rolled up graphene 
sheet, and interacting with EM-ficld in matrix represen- 
tation is 



The rotating wave approximation and the single-mode 
approximation of EM-field are used. The generalization 
for a multimode case thereafter being simple. All the 
components in §fi§ are considered to be hypercomplex 
operator n-numbers and they are the following. [Ha] 
represents the operator of cr-polaron subsystem energy in 
the absence of interaction between cr-polaron themselves 
and with EM-field. It is 



[# ff ] = (^o/2)£5>^.[ ei ] 



(7) 
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where a^ = \a m .., ..,,„, 



0"mj\ ~ \b m j) (b m j\ is z-transition op- 
erator between the ground and excited states of m-th 
cr-polaron in j'-th chain. The second term 
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[Hf] = fuo 2^ a + a[eiY 

3=0 



(8) 



[H] = [U a ] + [H F ] + [H aF ] + [H aa ] + [Hlf] 



(6) 



is the Hamiltonian of the free electromagnetic field, which 
is represented in the form of hypercomplex operator n- 
number. The component of the Hamiltonian © 

rs-l 

[U aF ] = hg Y, E «-« e4femQ + t-ja+e-^leiY 

(9) 
describes the interaction of cr-polaron sybsystem with 
EM-field, where o is the interaction constant. The Hamil- 
tonian 

n — 1 n— 1 

\u aa \ = -^E E 4-ii [ei]l E w ( fl ™+i^ m j 

1=0 j=0 m 

n — 1 n — 1 

- h zZ E 4-ii [ei]l E m (°»»-^i W 

1=0 j=0 rn 

n—1 n—1 

~ h zZ E 4-il N' E I W (Wi,l [ei] j 

1=0 j"=0 m 

n—1 n—1 

- h zZ E efwi n' E m fo-wi N J \ 

1=0 j=0 m 

(10) 

where ti&,2. \ are the energies, characterizing intrachain 
{I = j) and interchain {I ^ j) polaron-polaron interac- 
tion for the excited (£ < - 1 - > ) and ground (^ 2 - ) ) states of j'-th 
chain, \b m j) , \a m j) are ground and excited states corre- 
spondingly of m-th cr-polaron of j-th chain. Hamiltonian 
in the form like to (TTU)) at n = is usually used for de- 
scription of tunneling between the states with equal en- 
ergies, in particular for tunneling between quantum dot 
states [llj], m- Hamiltonian (ITU1) at any n describes ac- 
tually the connection between pairs of the states, which 
satisfy the following condition - the first state in any pair 
results from the second state (and vice versa) by time 
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Figure 1: Temporal dependence of the integral inversion in 
the QD chain for a coherent initial state of light ((n) — 4). 



Figure 2: Spectral dependence of the signal amplitude, corre- 
sponding to temporal dependence, presented in Fig.l 



reversal. It is known, that for given states Cooper ef- 
fect takes place, which, in particular, plays main role 
in BCS-theory of superconductivity. Therefore, the ap- 
plication of Hamiltonian like to (ITU1) is possible for any 
pair of time reversal symmetric states with equal energy. 
The component of the Hamiltonian [Hlf] corresponds 
to the local-field effects originated from the dipole-dipole 
electron-hole intrapolaron interaction and it is 



j—0 m 



(11) 

where P is dipole moment of cr-polaron, D is depolar- 
ization tensor. In (|7|) to (ITT1) [eip is j-th power of the 
circulant matrix [ei], which is 



[ei] 



1 


... 





1 ... 





... 1 


1 


... 



(12) 



By the way, if one omits last two terms in Hamiltonian 
©, it goes into n-chain generalization of well-known 
Tavis-Cummings Hamiltonian [3]. The solution with 
given hypercomplex Hamiltonian and the solution for one 
chain in [2[ will formally have the same mathematical 
form. Following to [2J, we will represent the state vector 
of the "NT+EM- field" system in terms of the eigenstates 
of isolated polaron and photon number states in the form 

[|*(*)>] = 

n-l 

T.iT.T.^ki^ i<w> +<*(*) i w>)}[eiF- 

j=0 l m 

(13) 



Here, \b mj ,l) = \b mj )®\l), \a mj ,l) = \a mj )®\l), where \l) 
is the EM-field Fock state with I photons, A 3 m , (£) , B 3 m l (£) 
are the unknown probability amplitudes. Then, solv- 
ing the nonstationary hypercomplex Schrodinger equa- 
tion for the state vector [|^(i))], we obtain in continuum 
limit 



[*W)] 



. n — 1 p 

[h lf ] = -£(?t?) J2 5>-^-> +^<^»n j ', J 



[W(h, 0)] exp{it([9 l (h)} - gVT+l[x])}e thx dh 



(14) 



where x is hypercomplex axis x — [x, x, ..., x], [<& (x,t)] is 

[^x^^exp ^-^K ^^ixM (15) 

in which [a z ] is Pauli z-matrix, A = r _1 is phenomenolog- 
ical relaxation factor. In its turn [\I/' (#,£)] is continuous 
limit of functional block matrix of discrete variable m, 
which is 



[*m,lW] 



[•Am,l(*)] 

[B m ,i+i{t)\ 



(16) 



that is, it is consisting of two [n x n] matrices of proba- 
bility amplitudes 



[a^(*)] = 5>U*)M J ', 

J=0 
n-l 

J=0 



(17) 



which are determined by relationship (|13[) . Further, ma- 



trix [9(h)] in flUJ} is 



n-l 



[8(h)] = ~{([0 1 (h)] + Mh)])®[E»]} 

+ l{(Mh)]-[e 2 (h)])®Wz]h 



(18) 



where [E%] is unit [2 x 2]-matrix, [9i(h)} and [02 (M] are 



[9i(h)] = [ti}{2-a 2 (h+~) 2 }, 

[9 2 {h)] = M{2-a\h-^f) 



(19) 



Here [£i], [£2] are [n x n] matrices of coefficients in ([TO]). 
that is 



ft — 1 ft — X 

$] = E #-j| N', Kal = E ft N', ( 2 °) 



;=o 



1=0 



where one takes into account, that in view of axial sym- 
metry [£j 2 ] do not depend on j. Consequently, we have 

[£{] = &]', [&] = [U Matrix [ X ] in © is 

[x] = [E n ] (g) [o-a,] exp-i[cr z ](w - w )t. (21) 

Matrix elements of [$'(£,£)] are 



& qp (x,t) = j 0g(/i,O)exp cxpihxx 



(22) 



exp{iNe x P O^'CO — gvT—lKj(h))}dh, 

3=0 



where 8'(/j,0),)Jj(/i), ^(/i) are determined by eigenval- 
ues of $'(/i, 0), 0(/i) and x(/i), which are considered to 
be n-numbers. They are 



1 1 

e'(M) = -k^a'CM)) = - y>J.(M)exp 



2nqji 



3=0 



1 



n-l 



0,-(fc) = -k,(0(h)) = - V r (fc) exp H-2I1. 
n n z -^ n 

r=0 



(23) 
(24) 



Kj(h) = -k,-(x(/0) = - E »W exp 



2njri 



(25) 



r=0 



Owing to axial symmetry, the hypercomplex solution 
obtained can be represented in the form of sum of n 
solutions for n chains, that is hypercomplex n-number 
&(x,t) is 



n-l 

9=0 

where the solution for q-th chain $' (x,t) is 



(26) 



-1 



.1 



p=0 



(x,t)[ei 



(27) 



in which the matrix elements $' „(x,t) are determined 
by (|22l) . The relationship (|25p by taking into account 
(fT4f - (l24l) determines Rabi-wave packet, which prop- 
agates along individual chain of zigzag NT. It is evi- 
dent, that the starting equation has also feasible solu- 
tion in the form of traveling Rabi-waves. It is also clear, 
that analysis of Rabi-wave packet dynamics for individ- 
ual NT-component will be the same, that in [2]. Let 
us reproduce Fig. 4 from |2J, where temporal dependence 
of the integral inversion in the QD chain for a coherent 
initial state of light ((n) = 4) is presented. It is seen, 
that an originally Gaussian packet temporally oscillates, 
at that oscillations collapse to zero quickly, but revive 
with time increasing in, and what is characteristic, in 
another area of space |2[- Temporal dependence (Fig.l) 
of the integral inversion gives in implicit form the way 
for comparison of theoretical results with any stationary 
optical experiments. Really, it is sufficient to make a 
Fourier transform of given temporal dependence (Fig. 2). 
It will be proportional to signal amplitudes of infrared 
(IR) absorbance, IR-transmittance, IR-reflectance or Ra- 
man scattering, since they are determined by population 
difference. Therefore, Rabi waves can be registered in 
ID-systems by conventional IR- and RS-measurements, if 
the electron-photon coupling is rather strong. The com- 
parison with experiment is given in [10( . 
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